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ABSTRACT 

The expected imaging capabilities of future Extremely Large Telescopes (ELTs) will offer the 
unique possibility to investigate the stellar population of distant galaxies from the photometry of 
the stars in very crowded fields. Using simulated images and photometric analysis we explore here 
two representative science cases aimed at recovering the characteristics of the stellar populations 
in the inner regions of distant galaxies. Specifically: case A) at the center of the disk of a giant 
spiral in the Centaurus Group, (/j,b ~ 21, distance of 4.6 Mpc); and, case B) at half of the 
effective radius of a giant elliptical in the Virgo Cluster (//i ~ 19.5, distance of 18 Mpc). We 
generate synthetic frames by distributing model stellar populations and adopting a representative 
instrumental set up, i.e. a 42 m Telescope operating close to the diffraction limit. The effect 
of crowding is discussed in detail showing how stars are measured preferentially brighter than 
they are as the confusion limit is approached. We find that (i) accurate photometry (p ~ 0.1, 
completeness > 90%) can be obtained for case B) down to I ~ 28.5, J ~ 27.5 allowing us to 
recover the stellar metallicity distribution in the inner regions of ellipticals in Virgo to within 
~ 0.1 dex; (ii) the same photometric accuracy holds for the science case A) down to J ~ 28.0, 
K ~ 27.0, enabling to reconstruct of the star formation history up to the Hubble time via simple 
star counts in diagnostic boxes. For this latter case we discuss the possibility of deriving more 
detailed information on the star formation history from the analysis of their Horizontal Branch 
stars. We show that the combined features of high sensitivity and angular resolution of ELTs 
may open a new era for our knowledge of the stellar content of galaxies of different morphological 
type up to the distance of the Virgo cluster. 

Subject headings: instrumentation: adaptive optics - stars: imaging - galaxies:stellar content 



1. Introduction 

One of the key issues in modern astronomy 
concerns the Star Formation History (SFH) in the 
Universe. Direct observations of galaxies up to 
high redshift can be used to map the SFH, but 
since the integrated galaxy light is dominated by 
the contribution of the most recent generations 
of stars, the information that can be derived on 
the underlying older stellar population is severely 



limited. However, the SFH in galaxies can be 
uncovered by interpreting the Colour-Magnitude 
Diagrams (CMD) of their stars, which contain 



the f o ssil record of the i r SFH (e.g.. iGreggio et al 



1 993 1 iHoltzman et all ll999L iBrown et all 1200 
Cole et al.ll2007ilMcQuinn et al.l2009llWeisz etal.1 



2011 : see also ICignoni fc TosH l20ld for a review 



This kind of studies require accurate photometry 
down to faint magnitudes in crowded fields, which, 
with current instrumentation, is feasible only for 
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the nearest galaxies. This implies a very limited 
sampling of the SFH in the Universe, with plenty 
of dwarfs, a few spirals and no giant elliptical. 

By the end of this decade this situation is ex- 
pected to change significantly, as a number of 
extremely large telesco pes, a s the Giant Magel- 
lan Telescope (GMT; I Johns! 120081). the T hirty 



Meter Telescope fTMT; ISzeto et all 120081 ) and 



the European Extremely Lar ge Telescope (E-ELT; 
Gilmozzi fc Spyromiliol 2007 ) could come into op- 
eration. The large collecting area of ELTs coupled 
with adaptive optics cameras able to deliver quasi- 
diffraction limited images will allow us to probe 
a wide volume, where we can access a significa- 
tive sample of galaxies. More importantly, it will 
be possible to study dense stellar fields. Indeed, 
because of crowding, stellar photometry in exter- 
nal galaxies is currently feasible only in regions of 
relatively low surface brightness. For giant galax- 
ies, this prevents us from deriving detailed SFH 
where most of the galaxy mass is. One of the key 
advantages of taking images with extremely large 
aperture telescopes is the exceptionally good im- 
age quality when the telescope can work close to 
its diffraction limit. This significant improvement 
has two fundamental effects : 1) it produces a dra- 
matic reduction of the background light over the 
point spread function (PSF) area considered for 
the photometry; 2) it allows a significant improve- 
ment in the spatial resolution (proportional to the 
telescope aperture). These combined advantages 
offer the unique opportunity to carry out observa- 
tions of faint targets in crowded and/or structured 
objects (galaxies) at large distance, that cannot be 
exploited by any other telescope of smaller aper- 
ture, neither ground nor space based. In summary, 
the excellent resolution capabilities of next gener- 
ation large aperture Telescopes will enable us to 
study high surface brightness regions, i.e. the in- 
ner parts of galaxies, where star formation was 
more conspicuous, as well as to address directly 
age and metallicity gradients. 

The photometric performance of ELTs oper- 
ating clos e to the diffraction limit wa s investi- 
gated by lOlsen. Blum fc Rigautl (|2003l ) and by 
Deep et al.l (|201ll ). Both studies focus on the im- 
pact of crowding conditions on the photometric 
quality, as quantified by the 1 a width of the er- 
ror distribution as a function of magnitude, for 
stellar fields with different surface brightness. In 



this approach errors are considered to be sym- 
metric, with stellar luminosities having the same 
probability of being overestimated or underesti- 
mated. However, crowding induces asymmetrical 
photometric errors, with an exce ss of stars mea- 
sured brighter than they are (e.g. 



Tosi et al.lll991l 



Gallart. Aparicio h Vilchezl Il996t iRenzinil Il99 



This results into artificially brightened features on 
the CMD which may induce a systematic error in 
those parameters (e.g. distances and stellar ages) 
which are derived from their luminosity. In addi- 
tion, crowding may also affect systematically the 
distribution of stars across the CMD, jeopardizing 
its interpretation in terms of star formation his- 
tory. The impact of this asymmetry is negligible 
under low-crowding conditions (e.g. when the 1<7 
width is smaller than about 0.1 mag), but grows 
rapidly with crowding. In this paper we present 
the results of end-to-end simulations of two spe- 
cific science cases aimed at investigating the SFH 
in galaxies, which fully exploit the unprecedented 
capabilities foreseen for ELTs. 

There are several interesting programs con- 
cerning the study of resolved stellar populations 
in external galaxies, like the analysis of multi- 
ple g enerations in Globular Clusters (see Piottol 
2009 and references therein) ; the determination of 



the init i al ma ss function and its variations (e.g. 



Kroupal 120011 ) ; the recovery of the spatially re- 



solve d SFH in isolated and in interacting galaxies 
(e.g. Cioni et"ahll2011 for the Magellanic system); 
tracing the galaxy assembly process through the 
counts of resolved star s and detection of streams 
(e.g. ( Ferguson! l2007h ). Here we focus on the 
possibility of deriving global information on the 
star formation history for a significative sample 
of galaxies with a modest time investment. This 
means targeting distant galaxies, so that the field 
of view samples a fair fraction of their mass, and 
observing galaxies in groups or clusters, so as to 
enable the study of a representative galaxy pop- 
ulation. On these premises we have concentrated 
on the following science cases: A) deriving basic 
information on the SFH in the central region of 
a disk galaxy at the distance of the Centaurus 
group; and B) studying the metallicity distribu- 
tion at half of the effective radius of an ellipti- 
cal galaxy in the Virgo Cluster. For these cases 
we have developed simulations of stellar photom- 
etry, based on synthetic stellar populations and 
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Fig. 1. — Expected limiting magnitudes for point 
sources observed with MICADO@E-ELT (solid 
lines) compared to the magnitude of RGB Tip, 
red HB and 10 Gyr old MS TO stars at a distance 
of 4.6 Mpc (dot-dashed lines). The limiting mag- 
nitudes, relative to a S/N of 20 (upper curve) and 
of 5 (lower curve), are computed for a total of 5 
hours integration time. 



assuming the expected performances of the Adap- 
tive Optics assisted n ear-IR camera (MICADO, 
Davies fe Genzelll2010l ) for the E-ELT (see details 
in Sect. 3.1). The resulting CMDs were the ana- 
lyzed to assess the impact of the photometric er- 
rors on the specific science goal. 

In section 2 we describe our selected science 
cases. Sect. 3 report detail of our simulation and 
the following photometric analysis. Results of this 
study are given in Sect. 4 with general conclusion 
summarized in Sect. 5. 



2. The science cases 

The derivation of the SFH in galaxies from 
the CMD of their resolved stars is nowadays 



method (Tosi et al. 


1991: Tolstov & Saha 


199a 


ADaricio et aljl996; Harris & Zaritskvl2001; 


Dolphin 


2002), which aims 


at minimizing the distance 



Fig. 2. — Synthetic CMD obtained with a con- 
stant rate of star formation over the last 12 Gyr, 
a Salpeter initial mass function (IMF), and a sim- 
ple age-metallicity relation which, starting from 
Z = 0, goes through the solar metallicity (Zq) at 
4.5 Gyr ago, and reaches a value of Z = 1.1 Zq at 
the current epoch. The simulation, computed with 
the YZVAR code by G.P. B ertelli with the 200 2 
version of the Padova tracks (jGirardi et al.ll2002h . 
contains 200000 stars brighter than Mk = —2, 
and corresponds to a total mass of formed stars of 
2.9 • 10 8 Mq. The four diagnostic boxes superim- 
posed contain 122 (RSG), 598 (BSG), 340 (AGB) 
and 22198 (RGB) synthetic stars. 



els by varying the strength of the successive 
star formation episodes. The method basically 
rests upon age-dating and counting sta rs in di- 
agno s tic regions on the CMD (see, e.g.. iGreggiol 



2002; iGreggio fe Renzi ni 2011), and the derived 



between the observed stellar density distribu- 
tion across the CMD and corresponding mod- 



SFH depends on the age sensitivity of charac- 
teristic features sampled in the o bserved CMD 
([Gallart. Zoccali fc Apariciol 120051 ). Stellar age 
dating is most accurate for stars at the main se- 
quence (MS) turn off (TO); at the other extreme, 
the red giant branch (RGB) collects stars born 
over (almost) the total range of ages present in 
the galaxy. Also the Horizontal Branch (HB) is 
populated by stars in a very wide range of ages, 
but the truly old populations can be distinguished 
as RR Lyrae stars or blue HB stars. 
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Fig. 3. — Age distributions of the synthetic stars 
in Fig. [5] which fall in the four diagnostic boxes, 
as labelled. Notice that each distribution is nor- 
malized to the total number of objects in the cor- 
responding box. 

Fig. [U shows the limiting magnitudes for 
isolated points sources that are expected for 
MICADO@E-ELT with a 5 hours of total inte- 
gration, together with the typical magnitudes of 
RGB Tip stars, red HB stars and MS TO stars 
with an age of 10 Gyr, at a distance modulus 
of DM = 28.3 Comparing the observational limits 
with these characteristic magnitudes we derive the 
typical distances up to which these features can 
be detected and measured with a relatively good 
S/N (i.e. ~ 10, see Table [T]). It appears that we 
could sample old MS TOs only within the Local 
Group and little beyond, while for galaxies in the 
Centaurus group (DM = 28.3) we will measure HB 
stars with high S/N. In the nearest galaxy cluster 
(Virgo, DM=31.3), the portion of CMD accessible 
with a reasonable exposure time is limited to the 
(upper) RGB and brighter magnitudes. Actually, 
HST already provides us CMDs down to these 
limits but only in sparsely populated regions; for 
example in Centauru s A, HB stars were detected 
Reikuba et al.1 (l2005h with ~ 9 hr exposures, in a 
region located at ~ 40 Kpc from the center (/zk ~ 
26). At this distance, the portion of the CMD at 
magnitudes brighter than Mi ~ — 2 has been mea- 



Table 1: Maximum distance moduli (DM) to ob- 
tain accurate J band photometry (S/N ~ 10) of 
isolated stars in various evolutionary stages, with 5 
hours of total integration with MICADO@E-ELT. 



Evolutionary 


Stage 


DM 


Old MS Turn 


Off 


26.5 


RR Lyrae 




27.5 


Red HB 




30.0 


tip of RGB 




34.5 



sured and analyzed for late typ e dwarfs around 
M83 (ICrnoievic. Grebel fc Colell201lh . Similarly, 
photometry down to ~ 1 mag nitude bel o w the Tip 
of the RGB was obtained by ICaldwelll (|2006l ) for 
dwarf galaxies in Virgo. The real gain we expect 
with MICADO concerns the ability of deriving 
data of similar quality in high surface brightness 
regions within galaxies, to the advantage of sam- 
pling the bulk of the star formation in medium 
size and giant galaxies. The simulations described 
in this paper aim at checking this quantitatively 
for the two specific science cases described below. 

2.1. The SFH in a disk galaxy in the Cen- 
taurus group 

Fig. [2] shows a synthetic CMD which can be 
considered as representative of the stellar popu- 
lation in the disk of a late type galaxy, with the 
synthetic stars color coded according to their age, 
as labelled. We examine here the possibility of 
recovering basic information on the SFH for such 
a stellar population placed at a distance of ~ 4.5 
Mpc. Four diagnostic boxes (see Table [2]) are su- 
perimposed on the synthetic CMD, with the aim 
of targeting specific age ranges, i.e. the Red Su- 
pergiant (RSG) box for the youngest component; 
the Blue Supergiants (BSG) and Asymptotic Gi- 
ant Branch (AGB) boxes for the intermediate age 
components; the RGB box for the old compo- 
nent. The age distributions of the stars falling 
in each of the four boxes is plotted on Fig. [3j 
showing that there is little age overlap in the var- 
ious boxes. The simulation allows us to estimate 
the average number of stars per unit stellar mass 
formed in the age range sampled by each box, here- 
after referred to as specific productions (Pbox), re- 
ported in Table [2] Star counts in these diagnostic 
boxes performed on an observational CMD allow 
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Fig. 4. — Synthetic CMDs fo r a stellar population with a flat age distribution between 10 and 12 Gyr, and a 
wide metallicity distribution (jReikuba et al.ll2005l ) plotted in the right panel. Dots are coloured according to 
their metallicity, with the same encoding as in the right panel. The simulation, computed with the YZVAR 
code by G.P. Bertelli with the 2002 Padova tracks database and adopting a Salpeter IMF, contains 200000 
stars brighter than Mr- = +1, and corresponds to a total mass of formed stars of 6.8 TO 7 M Q . The RGB 
diagnostic box on Fig. [5] is reproduced on the left panel. Notice that the color range of the infrared CMD 
(left panel) is much narrower than that of the optical- infrared CMD (central panel). 



us to recover quantitative information of the SFH 
in the observed field over the whole galaxy life- 
time: it suffice to divide the star counts by the 
appropriate specific productions to derive the av- 
erage mass turned into stars in the corresponding 
age range. This method is not as sophisticated 
as the full synthetic CMD method, which fits the 
detailed stellar distribution using small boxes, but 
it's much more straightforward and easily applica- 
ble to large datasets. At any rate, the reliability 
of the derived SFH depends on the completeness 
of the counts within the boxes, sensitive to the 
photometric error which moves the stars from or 
into adjacent regions of the CMD. Therefore we 
expect a dependence of the result on the crowding 
conditions, i.e. on the surface brightness of the 
sampled galaxy region. In the following we test 
how the star counts in these boxes are affected for 
a stellar field in the very center of a disk galaxy lo- 
cated in the Centaurus Group adopting a distance 
modulus on m — M = 28.3 and no extinction. 



Table 2: Properties of the diagnostic boxes for the 
YOUNG population science case: column (2) re- 
ports the age range sampled, column (3) the av- 
erage specific production (for a Salpeter IMF) in 
Mq , columns (4) and (5): number of synthetic 
stars in the input and output CMD of our simula- 
tion, which refers to a total stellar mass of ~ 2.6 
x 10 7 Mo in the FoV. 



Box 


Age Range 


^box 






Nout 


RSG 


< 40 Myr 


1.26 xlO" 


-4 


14 


14 


BSG 


( 40 to 120) Myr 


2.91 xl0~ 


-4 


59 


59 


AGB 


( 0.12 to 1) Gyr 


1.25 Xl0" 


-5 


35 


35 


RGB 


> 1 Gyr 


8.54 xlO" 


-5 


1933 


1929 



5 



2.2. The metallicity distribution for stars 
in a giant elliptical in Virgo 

Fig. @] shows the luminous portion of the CMD 
of an old stellar population with a wide metallicity 
distribution, akin to elliptical galaxies. Of the four 
diagnostic boxes drawn on Fig. [5] only the RGB 
box is populated, and includes 3232 stars, yield- 
ing a specific production of -Prgb = 4.75 10~ 5 Mq 1 
for this stellar population. The absence of stars in 
the RSG, BSG and AGB boxes implies (virtually) 
no star formation at ages younger than ~ 1.5 Gyr, 
whereas dividing the star counts in the RGB box 
by Prgb one obtains the mass transformed into 
stars in the field. More interestingly, the color 
distribution of the bright RGB stars can be used 
to recover their metalli city distribution, as in, e.g., 
Harris fc Harris! ([20021 ) . The different metallicity 
bins separate much better on the central panel of 
Fig. 3J due to the higher temperature sensitiv- 
ity of the wider color baseline. In this respect 
an optical-infrared color provides the best lever- 
age for the determination of the metallicity distri- 
bution of the RGB stars. The horizontal bars in 
the left and central panels of Fig. |4] show the size 
of a 0.1 mag error on the respective color: clearly 
the photometric accuracy needed to recover the 
metallicity distribution is less demanding in the 
optical-infrared CMD than in the infrared CMD. 
At any rate, photometric errors smear out the true 
color distribution of the stars, thereby affecting 
the derivation of the underlying metallicity distri- 
bution. Thus, the accuracy of this distribution 
determined from the analysis of the CMD is sen- 
sitive to crowding conditions. In the following we 
test the effect of the photometric errors for a stel- 
lar population located at half the effective radius 
of an elliptical galaxy in the Virgo Cluster, adopt- 
ing a distance modulus of m — M = 31.3 and no 
extinction. 

3. Simulations and Analysis 

Simulated images were constructed with a toofl 
developed during the Phase A stud y of MICADO 
(jFalomo. Fantinel fc Uslenghill201ll ). The simula- 
tor is able to generate synthetic frames with a large 
variety of telescope and instrument parameters as 




Fig. 5.— MAORY PSF in J band (FoV 0.05 arc- 
sec) MICADO@E-ELT. 



well as for input field of stars and/or galaxies. In 
the MICADO@E-ELT configuration all the rele- 
vant parameters of the telescope and of the in- 
strument are included, as well as the conditions of 
the observations. For the simulated images pre- 
sented here we assumed a 1277 m 2 collecting area, 
a total throughput of 0.4 and 0.39 respectively in 
the / and K bands, a read noise of 5 electrons, and 
a plate-scale of 3 mas per pixel. We assumed the 
typical near- IR sky background at Paranal (Chile) 
and included the contribution of thermal emission 
in the near- IR bands; specifically we adopted a 
background os 20.1, 16.3 and 12.8 mags in the /, 
J and K bands respectively. 

3.1. MICADO 

MICADO is the Multi-AO Imaging Camera for 
Deep Observations, designed to wor k with adap- 



tive o ptics (AO) on the 42m E-ELT ([Davies et al 



J The Advanced Exposure Time Calculator, AETC; 
http:/ /aetc. oapd.inaf.it/ 



2010l ). MICADO is optimized for imaging at the 
diffraction limit, and will fully sample the 6-10mas 
full width half maximum (FWHM) in the J-K 
bands. With a throughput exceeding 60% its sen- 
sitivity at 1-2 [im will be comparable to that of 
JWST (James Webb Space Telescope) for isolated 
point sources, but MICADO's exquisite resolution 
will allow us to investigate regions with a much 
higher crowding, compared to JWST. The instru- 
ment is able to image a relatively wide (53" x 53") 
Field of View, and a large number of wide and 
narrow band filters is foreseen, ranging from the / 
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RADIUS (arcsec) 



Fig. 6. — Radial growth of the encircled energy of 
the PSFs adopted for our synthetic frames. The 
inset shows the radii (in arcsec) which include 50 
% (R0.5) and 90 % (R0.9) of the total energy. 



to the K band. 

The simulations assume that the images are 
obtained using MICADO in conjunction with 
MAORY, the post-focal adaptive optics modules 
currently under study for the Europ ean Extremely 
Larg e Telescope. The M AORY (jDiolaiti et al. 
20101 : iFoppiani et al.l l201dT) module uses six laser 
guide stars to produce multi-conjugate AO (MCAO) 
correction. We assumed the PSF computed with 
a seeing of 0.6 arcsec d, with FWHM of 4 mas, 5 
mas and 9 mas in the /, J and K bands respec- 
tively. Since the Field of View of our simulation is 
small compared with the expected PSF variation 
on the total field of MICADO, we assumed a fixed 
PSF over the whole simulated frame. Specifically, 
we used the MAORY central PSF mapped over 
2048 x 2048 pixels. In Figure [5] we show a rep- 
resentative image of the adopted PSF in J band. 
The PSF in the I and K bands have a very sim- 
ilar shape but different Strehl (i.e. the ratio of 
peak diffraction intensities of an aberrated versus 
perfect wavefront) and encircled energy values. 
Fig. [6] shows the radial growth of the encircled 

2 http:/ /www. bo. astro.it/maory/Maory /Welcome. html 
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Fig. 7. — Central region (1.5 x 1.5 arcsec) of the 
simulated image of the disk galaxy (case A) in 
the J band as would be observed in 5 hr exposure 
using MICADO@E-ELT. The surface brightness is 
fij= 18.9 □". 



energy in the photometric bands used for our 
simulations. These distributions present a very 
sharp central peak, which becomes progressively 
more pronounced going from the / to the K band, 
reflecting the trend of the Strehl ratio. A quanti- 
tative characterization of the PSF is included in 

Fig. El 

3.2. Input stellar populations 

Table [3] lists the properties of the input stel- 
lar populations of the two science cases: A) the 
central region of a disk galaxy in the Centaurus 
group (hereafter YOUNG stellar population) and 
B) a region located at half of the effective radius of 
an Elliptical in the Virgo cluster (hereafter OLD 
stellar population). Illustration of CMDs for the 
two stellar populations are shown in Figs [2] and [H 
where we also specify the adopted star formation 
histories. 

The input stellar lists are generated from theo- 
retical simulations following a procedure described 
in detail in the Appendix; here we only mention 
that special care was devoted to ensure (i) ad- 
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Table 3: Parameters of the simulated stellar populations. Column (2): ratio between the total mass of 
formed stars and the -B-band light in solar units; columns (3) to (6): integrated colors; columns (7) to (9): 
adopted surface brightness, distance modulus and Field of View; column (10): total mass (in solar units) of 



the simulated stellar population following 


from the 


previous 


entries. 










Stellar Population M SF /L B B-V 


B- I 


B-J 


B-K 


Mb 


DM 


FoV 


M SF 


YOUNG 1.52 0.48 


1.29 


2.16 


2.99 


21.07 


28.3 


12" x 12" 


2.66 xl0 y 


OLD 7.05 0.88 


1.97 


2.84 


3.66 


21.6 


31.3 


3" x 3" 


7.53 xlO 7 





for* »*:v.^>-:^5&v..» • 



Fig. 8. — Central region (1.5 x 1.5 arcsec) of the 
simulated image of the elliptical galaxy (case B) in 
the J band as would be observed in 5 hr exposure 
using MICADO@E-ELT. The surface brightness is 
/ij = 18.8 □". 

equate sampling of the short lived evolutionary 
stages, and (ii) completeness of the input stellar 
list in all the photometric bands used in the spe- 
cific science case. The size of the simulated pop- 
ulation follows from the proportionality between 
the number of stars brighter than a given magni- 
tude limit and the total luminosity of the stellar 
population sampled by the synthetic frame: 



Lr = i7 ol /l -0-4*(MB-5.48-DM) r 



,0 



(1) 



where FoV is the Field of View in square arc sec- 
onds, ix-q is the (un reddened) surface brightness 
of the stellar population, DM is the distance mod- 



ulus, and we have adopted a solar absolute mag- 
nitude in the B band of Mb,© = 5.48. 

In order to keep the input stellar list within 
a manageable size only objects brighter than 
a threshold magnitude are considered as point 
sources, while the remaining light of the stel- 
lar population is distributed over the frame as 
a pedestal, with its associated Poisson noise. The 
threshold magnitude is chosen ~ 1.5 mag fainter 
than the limiting magnitude ( S/N=5 ) of the tele- 
scope + instrument combination for the assumed 
exposure time. This ensures that the effect of 
blending of stellar images is well described also at 
the faint end of the luminosity function. 

3.3. Synthetic Frames 

The AETC tool with the MICADO configura- 
tion was used to generate images in the /, J and 
K bands for the two science cases in Table |3] No- 
tice that the mass to light ratio and the colors are 
univocally determined by the adopted SFH, while 
the stellar mass (Msf)> sampled by the frame, de- 
pends on the chosen surface brightness, distance 
and FoV. For the YOUNG case, the input stellar 
list consists of 644483 stars brighter than K =30.7; 
for the OLD case, it counts 125058 stars brighter 
than _ftT=31. The input lists specify mass, age, 
metallicity, magnitudes in the /, J, and K bands, 
and coordinates on the frame for each star. Syn- 
thetic frames are generated from them in the cor- 
responding bands and Figures [7] and [H] show sec- 
tions of the J band simulated images. The wide 
magnitude range of the objects in Fig. [7] is read- 
ily visible, with a few very bright stars showing 
the characteristic pattern of the PSF. It can be 
noticed that many objects are detected under the 
halo of the bright stars due to the narrow core 
of the PSF. The image in Fig. [8] has a very dif- 
ferent appearance, with many stars of comparable 
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Fig. 9. — Input and output luminosity functions in 
the J (left) and K (right) bands for the YOUNG 
stellar population science case. The thin lines 
show the luminosity function of the whole input 
list. The shaded histogram highlights the luminos- 
ity function of those input stars which have been 
position matched in the reduction process. The 
thick lines show the output luminosity function. 
Notice the effect of blending particularly evident 
at the faint end, so that, e.g., stars with input 
magnitude as faint as K ~ 30 are detected in the 
reduction process, but are measured much brighter 
than they are. 



brightness: this reflects the different properties of 
the luminosity function of the two populations. 

3.4. Photometric Analysis and Methodol- 
ogy 

We performed the extraction of the P SF pho- 
tomet ry on the simulated frames using the [Stetson 
(|l994T) suite of programs: DAOPHOT, ALLSTAR, 
ALLFRAME, and related supporting software. 
The usual steps of detecting the sources in the 
images, creating the PSF and finally measuring 
the PSF photometry were executed via a semi- 
automated procedure. To avoid the introduction 
of any possible bias, the image simulations were 
performed by a person different from the one who 
carried on the photometry of the fields. There- 
fore the analysis was performed with no previ- 
ous knowledge of the simulated images. Only the 
"observational" parameters present in the header 
of the image has been used, like the number of 
summed and averaged exposures, the gain of the 
detectors, the read out noise etc. 

Performing PSF photometry with DAOPHOT 
requires two main steps: the construction of the 



PSF and the extraction of the photometric mea- 
surements. For the PSF creation we used 500 stars 
per image and we calculated all the 6 different op- 
tions for the analytical! modeling of the PSF, let- 
ting the program decide for the best PSF model 
to be used based on the overall minimum \ 2 ■ The 
procedure was iterated up to 5 times. At each it- 
eration we instructed DAOPHOT to subtract all 
the PSF neighboring stars from the original image 
and used this new image to extract the updated 
PSF model. We also discarded all PSF candidates 
with bright neighbors, or with saturated/dead pix- 
els. This procedure turned out very effective, typ- 
ically converging to a suitable PSF within 3 iter- 
ations. The particular shape of the MAORY PSF 
with its sharp core and extended wings forced us 
to use a DAOPHOT PSF dimension much larger 
than the one used for typical gaussian-like PSF. 
For the latter case the PSF radius of DAOPHOT 
is 4 -J- 5x the FWHM of the stars, while in the 
case of MAORY we had to use a factor of 15 x 
the FWHM of the core. In this way we enclose 
~ 95% of the stellar flux within the PSF radius. 
The analytic part of the of PSF was reduced to a 
minimum allowing for the DAOPHOT fitting ra- 
dius a value of ~ 0.9 x that of the core FWHM. 
The iterative procedure outlined above ensures an 
accurate description of the PSF, including its wide 
wings. The FWHM of the PSFs finally extracted 
by DAOPHOT are of ~ 1.8, 2.2 and 3.8 pixels 
in the /, J and K band images respectively, with 
a very small difference in the two science cases. 
These measured FWHM are only slightly larger 
than those of the input PSFs (see Sect. 3.1). Al- 
though we set a fixed (non position dependent) 
PSF, we allowed for internal variability of the re- 
constructed DAOPHOT PSF. 

The final extraction of the photometry is per- 
formed with the ALLSTAR+ALLFRAME pro- 
grams. ALLSTAR extract the PSF magnitudes 
on the image using the model PSF reconstructed 
by DAOPHOT. When more than one image of 
the same field is available (even in different fil- 
ters) the extracted photometry can be used as in- 
put to ALLFRAME, which repeats the ALLSTAR 
procedure simultaneously on all images. This con- 



3 DAOPHOT allows the calculation of 6 different analytical 
model PSF via its analytic function parameter: Gaus- 
sian, Moffat 1.5, Moffat 2.5, Lorentz, Penny 4 and Penny 
5 functions. 
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Fig. 10. — Photometric error as function of the input magnitude in the J (left) and K (right) bands of 
all matched sources for the YOUNG stellar population case. Dark (light) points are sources with a single 
(multiple) candidate counterpart in the input list. Solid lines show the la uncertainty for isolated stars. 
Dashed lines show the median error of the matched stars: notice that because of blending more than half 
of the stars are measured brighter than they are at J, K > 28. The filled circles show the la widths of the 
negative and of the positive error distributions. The median error and the +a width are plotted only up to 
the value of input magnitude for which the error distribution on the positive side is fully sampled. 



siderably improves the final photometric precision 
because of the much more significant statistical 
treatment of each detected star in the observed 
field, both position and photometry are being con- 
strained by all images. The use of ALLFRAME 
allowed us to extend the final list of the measured 
objects at a ~ 15 % fainter limit than the detec- 
tion threshold in a single band. 

The previously described procedure allowed the 
measurement of the PSF photometry of each im- 
age in an automatic, un-supervised mode requiring 
only one input parameter: the extraction thresh- 
old. We usually set this parameter to 3.0 x a of 
the background, while all other parameters were 
set internal to the procedure. The final product 
is a catalogue of magnitudes measured with ALL- 
FRAME operating on the J and K images for the 
YOUNG stellar population science case, on the /, 



J and K band images for the OLD one. 
4. Results 

4.1. Photometric Quality 

Having full knowledge of the input star list we 
can evaluate the quality of the photometry object 
by object. In order to pair detections to input 
sources it is necessary to define a matching algo- 
rithm, a task which is not completely trivial. We 
use a search radius of 1 pixel to select candidate 
mates of each detected source; within such area we 
may find (a) no input source, (b) one input source, 
(c) more than one input source. Case (a) corre- 
sponds to spurious detections, due to, e.g., noise 
spikes, or interference of the wings of the PSF of 
neighboring bright stars; case (b) corresponds to 
a one to one match, while in case (c) we have the 
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Fig. 11. — Input and output luminosity functions 
in the I (left) and K (right) bands for the OLD 
stellar population science case. Thin and thick 
lines show respectively the input and the output 
luminosity functions; the shaded histogram high- 
lights the luminosity function of those input stars 
which have been position matched. Notice that 
the large distance prevents us to recognize the 
populous HB/clump feature. 

embarrassment of more than one possible (input 
star) mate. We adopt the criterion of choosing as 
counterpart the brightest among the candidates. 
Different criteria could be envisaged (e.g. impos- 
ing a limit on the difference between the input 
and output magnitudes), all however getting into 
trouble at faint magnitudes, when approaching the 
confusion limit. Our choice has the merit of be- 
ing totally independent of the magnitude differ- 
ence between the input and output source, which 
is in fact the subject of the analysis described in 
the following. 

4-. 1.1. The Young Stellar Population at 4-6 Mpc 

The simulated images for this science case in- 
clude ~ 650000 stars with J,K< 31; the photo- 
metric measurement yields ~ 163000 sources de- 
tected in J and K of which ~ 159000 have a 
matched counterpart in the input list, while ~ 
4000 (i.e. ~ 2 % of the detections) are spurious 
objects. Most of these fake sources are located 
around input bright stars, and are features of the 
PSF (see Fig. |3J) misinterpreted as stars by the re- 
duction package because of a non perfect subtrac- 
tion of the PSF. Fig. [9] compares the luminosity 
function of the output catalogue (thick lines) to 
that of the input data (thin line). The luminos- 
ity function appears very well recovered brighter 



than J ~ 26.5, if ~ 25.5, while a strong incom- 
pleteness is evident fainter than J, K ~ 29. In 
the intermediate magnitude range the output lu- 
minosity function counts an excess of stars com- 
pared to the input data because of blending, i.e. 
a luminosity increase of the detected source due 
to overlapping PSFs, as argued below. In Fig. [9] 
we also show the luminosity function of the input 
stars which have been position matched (shaded 
histogram): clearly only a fraction of the input 
stars are recovered in each bin, which is the true 
completeness, or the probability of detecting a star 
of given input magnitude. Due to photometric er- 
rors, the position-matched stars can be found in 
different bins, typically brighter than their true 
value, which is the reason why the output luminos- 
ity function exceeds the input one at most mag- 
nitudes except for the brightest bins, where pho- 
tometry is very accurate. At the faintest magni- 
tudes incompleteness is boosted by the low S/N 
and stellar background, and the output luminos- 
ity function drops to zero. 

The photometric quality is further illustrated in 
Fig. [10] where we compare the magnitudes of the 
matched objects. The brightest stars are clearly 
recovered with a small error, but as the input 
magnitude increases, the distribution of the dis- 
crepancies widens as the stars are measured ei- 
ther fainter or brighter than they are. Close to 
the faint end of the input magnitude distribution 
the source detection is favored if the local back- 
ground is low, and/or the star is blended with 
other sources. Therefore, at the faintest input 
magnitudes, positive errors (i.e. output magni- 
tude ~> input magnitude) become under-sampled, 
and below some threshold, stars are recovered only 
if measured brighter than they are. The slanted 
limit in the right panel of Fig. [10] corresponds to 
the detection threshold used in the reduction of 
the K band image, i.e. an output magnitude of 
~ 29. The effect of blending is present over a 
wide magnitude range, and leads to an asymmet- 
rical distribution of the photometric errors, which 
extends over a wider range on the negative side 
compared to the positive side. Notice for exam- 
ple that the clump stars (at J ~ 27.3) can be 
recovered up to 1 mag brighter, or up to 0.5 mag 
fainter than they are. Although the average er- 
ror at this magnitude level is very close to zero, 
and the peak of the luminosity function due to 
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clump stars well recovered by the photometry (see 
Fig. E}, the spurious excess of stars above the 
true clump can be misinterpreted when analyzing 
the CMD, e.g. as due to a younger stellar gener- 
ation. Fig. [TU] also shows that negative errors are 
much more frequent for sources for which multiple 
candidate counterparts where found in the input 
stellar list, again demonstrating the constructive 
interference of stellar crowding, which leads to an 
artificial brightening of the detected sources. 

The asymmetry of the photometric errors is 
traced by the median of their distribution (which 
becomes progressively more negative as the in- 
put magnitude gets fainter), as well as by the la 
widths of the error distributions, computed inde- 
pendently for the stars measured brighter (lower 
locus in Figs [10] and [T2j and those measured 
fainter (upper locus) than they truly are. The 
asymmetry of the la width appears much more 
pronounced in the J rather than in the K band 
because of the higher background in the latter fil- 
ter. Actually the total error mainly results from 
the shot noise of the sky and thermal background 
plus the error associated to the underlying stellar 
population. In the K band the first contribution 
(which induces a symmetrical error) largely dom- 
inates over the second. The opposite holds in the 
J band, where the effect of the underlying stellar 
population can be better appreciated as an asym- 
metrical error distribution due to blending. 

4-1-2. The Old Stellar Population at 18 Mpc 

The simulated images for this science case in- 
clude ~ 125000 stars with K < 31 and J, I mag- 
nitudes down to 31.7 and 32.5 respectively. The 
photometric measurement yields ~ 22350 sources 
detected in /, J and K of which ~ 21700 have 
a matched counterpart in the input list, while ~ 
650 (i.e. ~ 3 % of the detections) are spurious 
objects. Also in this case, most of the spurious 
detections come from a non perfect PSF subtrac- 
tion, and also in this case the spurious sources are 
a small fraction of the total number of objects, so 
that the overall appearance of the CMD and of the 
LF is not significantly affected. 

As for the YOUNG stellar population, Figs HQ 
and [12] illustrate the photometric quality for the 
OLD stellar population. This science case differs 
from the previous one in the distance modulus (3 
magnitudes fainter) and in the magnitude range 
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Fig. 13 — Input (a) and output (b) CMDs for 
the YOUNG stellar population science case. The 
four diagnostic boxes are the same as in Fig. [2] 
The J = 26.5, 28 (dashed lines) and K = 28 (solid 
lines) loci are drawn to better appreciate the effect 
of the photometric error on the CMD features. 

covered by the input star list, because of the ab- 
sence of the young and intermediate age compo- 
nents. The input luminosity function (thin lines 
in Fig. HU does include the HB, at 31 < I < 32, 
29.5 < K < 30.5, but only very few of its stars 
are position matched (see shaded histogram), and 
in most cases they are assigned a too bright mag- 
nitude. The migration of stars towards brighter 
bins along the luminosity function due to blend- 
ing is evident also in this case, and can be appre- 
ciated as the difference between the shaded region 
and the output luminosity function (thick line his- 
togram). It is interesting to note that the cut off 
of the I band output luminosity function is more 
smooth than that of the K band one. This is 
due to the measurement procedure adopted which, 
once a source is detected on the K image, forces 
the detection at the shorter wavelengths, so that 
the adopted detection threshold translates into a 
sharp cut off only in the K band luminosity func- 
tion. 

The distribution of the photometric errors (Fig. 
1X2"]) presents similar characteristics to those de- 
scribed for the YOUNG stellar population case: 
due to blending, the median error is negative al- 
most all along the magnitude range covered, and 
the error distributions are skewed towards the neg- 
ative side. Compared to the YOUNG one, the 
OLD population science case is characterized by a 
much higher crowding, in spite of the similar sur- 
face brightness, because of the ~ 4 times larger 
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Fig. 12. — Photometric errors as function of the input / (left) and K (right) magnitude, for sources with 
single match (dark points) and with multiple match (light points) on the image. The dashed line shows the 
median error; the big dots show the ler widths of the error distributions computed separately for positive 
and negative errors. Notice that clump stars with input magnitudes / ~ 31, K ~ 30 are identified only if 
their luminosity is boosted by blending, i.e. only with large and negative errors. 



distance. For this reason, the error distributions 
are wider for the OLD than for the YOUNG case, 
particularly on the negative side (see Figs [10] and 
Q21 right panels). 

The / band appears to be heavily affected by 
crowding (see Fig. [T2] left panel) : in spite of the 
low sky background in the optical, the ler loci 
show that the total error distributions are very 
wide, and skewed on the negative side, again show- 
ing the effect of stellar blending. We notice that 
the brightest sources are affected by a large er- 
ror, which exceeds what is expected from the shot 
noise. Rather than to the underlying stellar pop- 
ulation, we attribute this effect to a non optimal 
sampling of the PSF in the I band, which has a 
nominal FWHM of less than 2 pixels. Although 
DAOPHOT is able to reconstruct such peaked 
PSF, the measurement of the magnitude is compli- 
cated by uncertainties in the position of the cen- 
troid. The use of ALLFRAME, i.e. taking ad- 



vantage of the information on the J and K band 
images, greatly improves the quality of the / band 
photometry; nevertheless, an error associated to 
the under-sampling is present all over the magni- 
tude range, and particularly visible at the bright 
end. 

A science case similar to the one discussed 

albeit with 
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here is presented in lDeep et al. 
some differences, including the CMD of the syn- 
thetic ste llar populat i on, an d the exposure time, 
for whi ch iDeep et al l (l201lh adopt 1 hr. In addi- 
tion, in Deep et al.l (|201ll Vs study, the simulated 
PSF has a coarser sampling than ours, and the 
instrument background is neglected. In spite of 
these differences, for the same surface brightness 
our results are broadly compatible with theirs in 
the / and J bands, with similar 1 a values at the 
bright end. 

To summarize: the comparison of the input and 
output catalogues shows that, roughly for both 
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Fig. 14.— Input (left) and output (right) CMDs 
in the region of the HB/Clump for the YOUNG 
stellar population science case. The color encodes 
the star's ages, as labelled. 



science cases, the output luminosity functions be- 
come 50 % incomplete at J,K ~ 29, I ~ 30.5, 
while a photometric error of a ~ 0.2 mag is 
reached at J ~ 29.5, K ~ 28 and / ~ 30, sim- 
ilar to the S/N=5 limits in Fig. Q] except for the / 
band where the error is much larger due to crowd- 
ing 0. These are however only indicative figures 
because the error distribution is very asymmet- 
rical due to blending, an error of a ~ —0.2 mag 
being reached at J ~ 28.5, K ~ 27 and I ~ 29. As 
a result, the output luminosity function is affected 
by an overall migration of objects from fainter to 
brighter bins, and in order to assess the feasibil- 
ity of a science case it is necessary to perform the 
whole test, proceeding, in our case, with the anal- 
ysis of the output CMDs. 

4.2. Comparison of CMDs 

4-2.1. Disk galaxy at 4-6 Mpc 

Fig. [13] shows the input (panel (a)) and the 
output (panel (b)) CMDs for the YOUNG stellar 
population science case, with the diagnostic boxes 
superimposed. The characteristic plumes of the 
input CMD are well recognizable in the bright por- 
tion of the output diagram, but become progres- 
sively blurred towards fainter magnitudes. Three 
lines are drawn on the CMDs to distinguish three 
regimes: brighter than J ~ 26.5 the main fea- 
tures of the CMD in panel (a) are well reproduced 
in panel (b), including their narrow color width. 
Fainter than K ~ 28 (slanted line) the star's dis- 



4 Notice that J A B ^ J + 0.9, K AB ~ A" + 1.9, Iab ~ / + 0.1 



tribution on panel (b) is completely dominated by 
the photometric errors, and presents no resem- 
blance to the corresponding region in panel (a). 
In the range 26.5 < J < 28 the CMD in panel (b) 
still retains the main features of the correspond- 
ing region in panel (a), but the color distribution 
of the stars is wider, showing that the photometric 
errors heavily affect the appearance of the CMD. 
Therefore, in order to derive the SFH from the 
analysis of this region of the CMD one needs to 
quantitatively account for the photometric errors. 
Notice that this region includes the clump stars. 

Table [2] reports the star counts in the four di- 
agnostic boxes in the input and output CMDs, 
showing that the counts are unaffected by the pho- 
tometric error. Had we considered a more dis- 
tant galaxy, or a stellar field with higher surface 
brightness, the diagnostic boxes would be closer to 
the confusion limit and the stellar blending would 
affect significantly the co unts (starting from the 
RGB box, cf. Eq. (8) in lOlsen. Blum fc Rigaut 
20031) . This would result into a systematic error 
in the derived star formation. For the science case 
considered here, i.e. the central part of a disk 
galaxy at a distance of 4.6 Mpc, the diagnostic 
boxes are bright enough that the mass transformed 
into stars in the relative age ranges can be robustly 
derived. Notice that the whole FoV of MICADO 
will be almost 20 times wider than simulated here, 
providing many stars for a robust statistics. 

By applying the full synthetic CMD method to 
the simulated data in Fig. [13] a more detailed SFH 
can be derived, in particular disentangling the con- 
tribution of the various age components in the 
range ~ 1 to 12 Gyr. To this end, the HB/clump 
region may provide an effective diagnostic. The 
left panel of Fig. [14] shows the input CMD in the 
clump region with stars color coded according to 
their age bin. The different colors appear sepa- 
rated on the input CMD, where regions occupied 
by stars in different age bins can be distinguished: 
only objects older than 10 Gyr are found in the 
blue part of the horizontal branch; stars with age 
between 1 and 2 Gyrs populate the faintest enve- 
lope of the clump, and the typical magnitude and 
color of core helium burning stars is a function of 
the age. Notice that the various age components 
partially overlap, and it is not possible to con- 
struct clean diagnostic boxes similar to those used 
in the brightest part of the CMD. In addition, pho- 
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tometric errors substantially smear out the color 
distribution of the stars. The right panel of Fig. 
[Ftlshows the same section of the output CMD with 
the detected objects color coded according to their 
age. To build this diagram we have assigned to 
each output source the age of its input star mate 
according to the position matching criterion. The 
smearing effect of the photometric error is read- 
ily apparent: the different age components of the 
stellar population are very hard to distinguish, so 
that the result of the application of the full syn- 
thetic CMD method will be very sensitive to the 
age binning and to the photometric errors. 

J^.2.2. Elliptical galaxy at Virgo 

Fig. [T5l compares the input and output CMDs 
for the OLD stellar population science case. The 
number of stars falling in the RGB box is 3538 in 
panel (a) to be compared to the 3616 star counts 
in the same region of the CMD in panel (b). Due 
to blending, some stars just fainter than the lower 
limit of the box become brighter and shift inside 
the box, contaminating the counts. However, the 
effect is quantitatively very small, and the total 
mass in stars can be recovered using the theoreti- 
cal specific production within few percent. Over- 
all, the CMD of the detected sources appears very 
similar to the input CMD, especially brighter than 
J ~ 28. Our aim is to check the impact of the pho- 
tometric errors on a widely used method to derive 
the metallicity distribution of an old stellar popu- 
lation, which is based on the analysis of the CMD 
of bright RGB stars. 

Fig. [15] shows the upper RGB of the input 
(panel (a)) and output (panel (b)) stellar lists, 
with the color encoding the metallicity bin. The 
[Fe/H] of each detected source is assigned equal to 
that of its input star mate on the image. The dif- 
ferent metallicity bins appear well distinguishable 
on panel (b), especially in the upper 1.5 mag of 
the RGB. The derivation of the metallicity distri- 
bution from the analysis of the CMD is subject to 
uncertainties related to the age-metallicity degen- 
eracy and to the AGB contributio n to the counts 



in th is part of the CMD (see, e.g.. iReikuba et al 
2010). Besides these systematic effects, the photo- 
metric errors introduce an additional uncertainty, 
which we evaluate here. We have selected from 
the output catalogue all the stars brighter than 
J = 28, and assigned to each of them a metallic- 
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Fig. 15.— Input (a) and output (b) CMDs for the 
OLD stellar population science case. The color 
reflects the metallicity bin of the object with the 
same encoding as on Fig. |4] The metallicity of 
the output stars is identified on the basis of the 
positional coincidence with input objects on the J 
band image. The J = 28 line is drawn to better 
appreciate the effect of the photometric errors on 
the width of the RGB. 




Fig. 16. — Difference between the metallicity de- 
rived from the position of the detected source on 
the CMD and that derived from the position of 
the detected source on the frame. The upper panel 
shows the results from the analysis of the ( J, / — J) 
CMD; the lower panel from that on the (J,J — K) 
CMD. 
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ity equal to that of the star in the input catalogue 
which has the closest position on the CMD. The 
difference between this value of [Fe/H] and that 
assigned from the source position on the image 
can be viewed as the error on the metallicity de- 
rived from the analysis of the CMD due to mere 
photometric errors. Fig. IT61 shows this difference 
as a function of the output magnitude in the J 
band, the two panels plotting the results from two 
different color combinations. The photometric er- 
rors introduce a mild uncertainty on the determi- 
nation of [Fe/H] when using the (J, I - J) CMD: 
for most points the difference between the photo- 
metrically determined and the true metallicity is 
smaller than 0.1 dex. The error is larger when us- 
ing the (J, J — K) CMD, as anticipated in Sect. 
2.2. Fig. IT7l illustrates the the impact of the errors 
on the photometrically derived metallicity distri- 
bution, The two histograms appear very similar: 
both peak at [Fe/H] ~ -0.2 and have similar widths 
and shape. We notice that the photometrically de- 
rived distribution is slightly overpopulated on the 
low metallicity side of the peak. Overall, however, 
the two distributions are very close to each other, 
and we conclude that this science case is well fea- 
sible in inner regions of giant elliptical galaxies at 
the distance of the Virgo Cluster. 

5. Summary and Conclusions 

In this paper we have explored the feasibility 
of two specific scientific applications for the study 
of the SFH in distant galaxies with the expected 
performance of the 42 m E-ELT equipped with the 
MICADO camera working close to the diffraction 
limit of the telescope. We have focussed on gi- 
ant galaxies located in the Centaurus group and 
in the Virgo Cluster. These are regions of the 
nearby Universe where a variety of galaxy types 
and morphology can be found, so that a significa- 
tive sampling of the SFH in the Universe can be 
gathered. In spite of the large collecting area of 
E-ELT, old MS turn-offs cannot be investigated 
at these distances; however interesting informa- 
tion on the SFH can be derived from the analysis 
of the luminous portion of the CMD. The crucial 
advantage offered by the E-ELT working close to 
the diffraction limit is related to the exquisite spa- 
tial resolution that enables accurate photometry in 
crowded stellar fields. This will allow us to study 
the SFH in high surface brightness regions of giant 
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Fig. 17. — Distribution of [Fe/H] as derived from 
the positional coincidence on the J band frame 
compared to that derived from the nearest neigh- 
bor match on the (J, I - J) CMD. 

galaxies where most of the stellar mass is located, 
and to address directly stellar population gradi- 
ents, from the outskirts down to the inner regions 
of galaxies. 

We have investigated how photometric errors 
impact on two science cases: A) the study of the 
SFH at the center of the disk of a spiral galaxy at 
a distance of 4.6 Mpc (/xb — 21); B) the recovery 
of the metallicity distribution in the inner regions 
(at 1/2 of the effective radius) in an giant ellipti- 
cal galaxy at a distance of 18 Mpc (fi\ ~ 19.6). To 
this end we have produced synthetic frames in the 
/, J and K bands by distributing theoretical stel- 
lar populations constructed with suitable star for- 
mation histories, and adopting the PSF currently 
foreseen for the MICADO camera assisted by the 
MAORY adaptive optics module. The synthetic 
frames were measured with a standard package 
for crowded fields photometry to yield output cat- 
alogues of detected stars. The input and output 
stellar lists have been compared to determine the 
photometric quality, and the output CMDs have 
been analyzed to assess to which extent the scien- 
tific aims of the two cases can be realized. 

For both science cases, and in all photomet- 
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ric bands, blending of stellar sources leads to an 
asymmetrical error distribution, and a general mi- 
gration of star counts along the luminosity func- 
tion towards the brighter bins. Over a wide mag- 
nitude range we detect an excess of stars recov- 
ered brighter than what they truly are, and the 
distributions of the (output-input) magnitude dif- 
ference is more spread out on the negative (with 
respect to the positive) side. This effect becomes 
particularly relevant as one approaches the con- 
fusion limit, where stars are measured only if 
their luminosity is artifici ally boosted by blending 
(Gregeio k RenzinillioiH ). At the same time, as 
the confusion limit is approached, a progressively 
larger incompleteness affects the luminosity func- 
tion, partly balancing the excess due to blending. 
As a consequence, the measured luminosity func- 
tion does not appear very different from the input 
one (see Figsl9landfTT|). Nevertheless the asymme- 
try of the error distribution persists and can lead 
to misinterpretations of the CMD, e.g. distances 
from RGB tip (and HB) stars may be underesti- 
mated; too young ages may be inferred from Main 
Sequence turn off stars. Similarly, the age distri- 
bution from the core Helium burning/clump stars 
could result too skewed at the young end, and 
bright RGB stars overshooting the tip could be 
mistaken for relatively young AGB stars. Besides 
the systematic shift to brighter magnitudes of im- 
portant diagnostic features, the asymmetry of the 
error distribution affects in a systematic way the 
star counts in the boxes used to derive the SFH. 
This effect, which is very important when trying 
to derive information from stars with a magnitude 
close to the confusion limit, may also hamper the 
results in less severe crowding conditions. Here 
we have checked that for our cases A) and B), the 
photometric accuracy is sufficient for their scien- 
tific aims. 

Specifically, our results show that in the very 
central parts of disks of giant spirals located in 
the Centaurus group (case A): 

• accurate photometry (a ~ 0.1) can be ob- 
tained down to J - 28.0 and K ~ 27.0. At 
these magnitudes the luminosity functions 
are > 90 % complete; 

• the characteristic plumes of the input CMD 
are very well reproduced in the output CMD 
down to J ~ 26.5, and the star counts in 



bright diagnostic boxes which target specific 
age ranges are perfectly reproduced. This 
implies that the star formation history in 
the last ~ 1 Gyr can be accurately derived, 
and that a robust estimate of the mass trans- 
formed into stars at ages older than ~ 1 Gyr 
can be obtained from the star counts on the 
upper RGB; 

• the analysis of the clump/HB stars to re- 
cover the SFH prior to ~ 1 Gyr ago with a 
better age resolution requires the application 
of the full synthetic CMD method, which in- 
cludes a careful modelling of the photomet- 
ric errors. For the surface brightness consid- 
ered, the smearing of the distribution of the 
clump stars is significant, hampering a ro- 
bust solution. This case will be better per- 
formed in less crowded regions of disks. 

The simulations of the case B) stellar popula- 
tion shows that for giant ellipticals in the Virgo 
cluster: 

• accurate photometry (a ~ 0.1) can be ob- 
tained down to I ~ 28.5 and J ~ 27.5 at half 
of the effective radius (/xb = 21.6). At these 
magnitudes the luminosity functions are > 
90 % complete; 

• the observed CMD well reproduces the 
width of the RGB in the upper 2 magni- 
tudes, and the color distribution of bright 
RGB stars is only mildly altered by the pho- 
tometric errors, so that its interpretation in 
terms of metallicity distribution of the stel- 
lar population is only slightly affected. We 
estimate that the photometric errors intro- 
duce an uncertainty of ~ 0.1 dex on [Fe/H] 
when using the I — J color combination as 
a tracer for the effective temperature. 

Given the good photometric quality in the 
bright RGB region for science case B), we argue 
that it will be possible to study the SFH in the 
central parts of disks in spiral galaxies in Virgo 
from star counts in bright boxes, as those shown 
on Fig. This will allow us to contrast the 
old (prior to ~ 1 Gyr) to the more recent star 
formation activity over entire galaxies and across 
the Hubble sequence for all cluster members. The 
good photometric quality in the upper RGB will 
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also allow us to derive the metallicity distribu- 
tion almost all over giant galaxies, thereby tracing 
metallicity and population gradients. 

For disks galaxies in the Centaurus group, the 
analysis of the Red Clump/HB feature will pro- 
vide a more detailed information on old star for- 
mation episodes. The central regions of disks ap- 
pear too crowded for this purpose, but the results 
presented here have been obtained with a stan- 
dard package for the photometry, which was not 
designed to work with such a structured PSF as 
the one provided by MAORY. Likely, the photo- 
metric accuracy will improve when using reduc- 
tion packages more suitable for complex PSF that 
is produced by adaptive optics systems. 

Further simulations that include a variable PSF 
shape over the whole field of view and other effects 
(non homogeneous background, field distortions, 
etc.) are required in order to fully characterize 
the capabilities of ELT imaging under different 
conditions. At the same time, simulations con- 
structed with different choices for the stellar pop- 
ulations and crowding conditions are necessary to 
adequately explore the scientific return. These is- 
sues will be the subject of a forthcoming paper. 
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A. Procedure to generate the synthetic stars 

The input lists of the synthetic stars were generated by random extractions of individual objects which 
belong to a library of theoretical simulations previously computed. For each of the two kinds of populations 
in Table El four simulations were constructed (G.P. Bertelli, priv. comm.), each constrained to have 200000 
objects brighter than a given magnitude Afx,max = +5, +1, —2 and -4. The brighter the magnitude limit, 
the more massive the simulated stellar population, and shorter lived evolutionary phases become better 
sampled. Figs [5] and 2] show one of these mother simulations respectively for the YO UNG and for the OLD 
population. 

The four mother simulations allow us to construct an input stellar list for YOUNG and/or OLD stellar 
populations of any prescribed size, since the number of stars in any region of the CMD is proportional 
to the total light (or total mass) of its parent stellar population. In particular, the number of stars with 
magnitude in a given range Mk,x < Mr < Mk,2 is proportional to the total luminosity sampled. Therefore, 
once the size of the stellar population to be generated is given, e.g. its Lb through Eq. (TTJ), the number 
of objects on its CMD with Mk.i < Mk < Mk,i is determined by scaling from the mother simulation 
of adequate depth, and the procedure is repeated to cover the whole magnitude range. The input stellar 
list is thus constructed by composing partial CMDs, each containing objects randomly extracted from the 
mother simulations. Notice that the partial CMDs are limited in the if-band because so are the original 
mother simulations. In the input stellar list each object is a synthetic star, i.e. characterized by mass, age, 
metallicity, and magnitude in all photometric bands. 

In principle, the input list could contain stars down to Mk = +5, which is the limit of the deepest mother 
simulation. In practice, there is no need to generate a list down to the faintest magnitudes, which, unless the 
sampled luminosity is very low, would consist of an unmanageable large number of entries. Indeed, the stars 
on the synthetic frame will be measurable only down to some limiting magnitude function of the exposure 
time (for the given observational set-up). Therefore the stellar population on the frame is split into two 
components, i.e. individual stars brighter than Mx.iim, plus a pedestal of stellar light made of the stars 
fainter than this limit. The latter component is evenly distributed over the frame, with its Poisson noise. 

Special care has been devoted to ensure completeness of the list of synthetic stars in all photometric bands, 
since this list is the input to generate frames at different wavelengths. The extraction procedure described 
above necessarily generates a list which is complete in the K band with the risk of under-sampling the blue 
stars. Therefore the value of MK,iim is derived imposing completeness in the bluest band of the CMD to be 
constructed. In addition, the limiting magnitude for the input stellar list is chosen ~1.5 mag fainter than 
the S/N=5 level, to account for the effect on photometry due to stellar blending which leads to an artificial 
brightening of the sources. 

We are indebted to G.P. Bertelli for computing the synthetic stellar populations used in our simulations. 
Preliminary results of this work are included in the Scientific Analysis Report document of the Phase A 
study of the MICADO camera. We thank R. Bedin and E. Held for their help in an early stage of this 
investigation, and A. Renzini for many useful discussions. This work was supported by funding of the phase 
A study of MICADO, and by INAF TECNO funding AO@SW 1.05.01.17.06. 
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